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Abstract 

A general, user oriented coapirter orograa, 
called VNAP2, has been developed to calculate high 
Reynolds nuaber, internal /external flows. VNAP2 
solves the two-dlaenslona 1, tlae-dependcnt Navler- 
Stokes equations. The turbulenct. Is aodcled trlth 
either a alxlng-length, a one transport equation, or 
a two transport equation aodel. Interior grid 
points are cowputed using the explicit MacCoraack 
scheae with special procedures to speed up the cal- 
culation If the floe grid. All boundary conditions 
are calculated using a reference plane chara«.terls- 
tlc scheae with the vlscoua terms treated as source 
terms. Several Internal, external, and Intemal/ex- 
temal flow calculations are presented. 

In troduction 

The coaputatlon of high Reynolds nuaber flows 
has become a major tool In the analysis and design 
of aerospace vehicles. While Navler-Stokes solu- 
tions for coaplete vehicle configurations are still 
beyond the limits of present-day computers, corputa- 
tlonal techniques are used routinely In the analysis 
and design of various Individual components, e.g., 
airfoils, wing-body coablnatlous , Inlets and nox- 
tles. Host of these analyses use either purely In- 
vlscld or the so-called patched vlscous-lnvlscld 
techniques due to their greater coaputr.tlonal effi- 
ciency and ease of use over the more exact Navler- 
3tok.es solution methods. These approximate tech- 
niques often yield results of surt>rlsingly high ac- 
curacy.^ However, their use Is generally limited to 
problems involving wesk vlscous-lnvlscld intertc- 
tlons or to strong Interactions that are sufficient- 
ly well understood to be modeled empirically. 

Navler-Stokes (N-S) solution methods, on the 
other hand, sre not subject to these fui.<Iamental 
limitations and are applicable to more general 
classes of flow prcLleas. However, Che N-S methods 
haw ihot found widespread use for design purposes 
due primarily to their expense and difficulty of 
use. Computer tun times of several hocrs a^s not 
uncomiaon In so' /Ing a hlgh-Re>nolds number problem 
In which the viscous layer must be wrll reaolyed. 
Furthermr^re, application of N-S methods Is often 
viewed by the engineer as an arc requiring extensive 
knowledge of the nuaerlcal algorithm and considera- 
ble trial and error tj obtain a correct, converged 
solution. The latter Is often the result of at- 
tempting to use a N-S computer code which has been 
written to solve a very specific class of problems 
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and any not be sufficiently general, well-docuaent- 
ed, and user-orlenced. Clearly, as more efficient 
N-S algorlthiM and larger, faater coaputera become 
available, more attention «ist be given to cne de- 
velopment of user-oriented N-S codes If they are to 
receive prictlcml application. 

The purpose of this paper Is to describe a N-S 

computer program, VNAP2,^ which has evolved over a 
period of several years for solving a relatively 
v’ido class of steady and unsteady. Internal and ex- 
ternal flow problems. VNAP2 Is a acdlflad version 

of VNAP^ and solves Che cwo-dlaenalonal (axlsvmmsc- 
rlc), Clae-dependent, compressible Navler-Srokes 
equations. Both sli.gle and dual flowing streara nay 
be solved. The flow boundaries may be arbitrary 
curved solid walxs. Inflow/outflow boundaries, or 
free- jet envelopes. Turbulent as well as laminar 
and lavlscld flows any be treated. Some typical In- 
ternal and exterrwl flow geometries that may be 
solved are shown In Fig. 1. Although the VNAF2 code 
has been applied mainly to oosxle and Inlet flews, 
the relatively general treatment of geometries and 
flow boundaries allows a variety of other problems 
to be solved, e.g., airfoils, flow-through nacelles, 
free-shear flown and free- jet expansions. 

In this paper, the aethedoiogy used In develop- 
ing VNAF2 as a user-oriented, production-type compu- 
ter program la presented and acme results obtained 
In solving a variety of flow problems are shown. 
Although the results. In some cases, point out the 
need for further laprovcaents in both nuaerlcal and 
physical modeling, they do Illustrate, tc a great 
extent, that practical Navler-Stokes applications 
are possible. 

Cc vemlng Equations 

The VNAP? code solves the two-dlaenslonal (ax- 
Isymmetric), time-dependent, Navler-Stokes equa- 
tions. The turbulence Is modeled using either \ 
elxlng-length, a one transport equation or a two 
tranaport equation aodel. The alxlng-length aodel 

employs the Leunder et al^ model for free shear lay- 
ers end the Cebecl-Salth^ aodel for boundary layers. 
The one equation aodel Is that, of Daly^ while the 

two aquation aodel Is the Jones -Launder^ aodel. 
For details of these turbulence models. Including 
boundary conditions, see Ref. 2. VNAP2 employs an 
explicit artificial vlscoalty to stabilize the cal- 
culations for knock waves. This artificial viscosi- 
ty Is used In place of the fourth-order smoothing 

usually employed by MacCoraack . ^ ^ For details of the 
governing equations see Ref. 2. 
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Phy»lc«l ano Coapuf tlonal Flow Sp«c<». 

The phyalcal flow space ^eoaecry Is shown In 
Fig. 2. The tlow Is froa left to right. The upper 
boundary, called Che wall, can be either a solid 
boundary, a free Jet boundary, or an arbitrary sub- 
sonic (not sal CO the boundary) Inf low/oucf lew bound- 
ary. The lower boundary, called the center body, can 
be either a solid boundary or a plare (line) of sy*- 
aetry. The geoaecry can be either a single flowing 
streaa or. If the dual flow space >^lla are present, 
a dual flowing streaa. Tne dual flow space walla, 
shown In Fig. 2, aay begin In Che Interior and con- 
tinue to Che exit (Inlet geoaecry), aay begin at Che 
Inlet and teralnate In Che Interior (afterbody geoa- 
etry), as shown lit Fig. 2, or aay begin and end In 
Che Interior (airfoil geoaecry). All of the above 
boundaries aay be arbitrary curved boundaries pro- 
vided the y coordinate Is a single value function of 
X. Th.s single value function of x requires dual 
flow space walls, that begin or teralnate In the In- 
terior, do so with pointed ends. The points can be 
very blunt, but cannot be vertical walls The left 
boundary Is a subsonic, superson'c, or alxed Inflow 
boundary while the right boundary Is a subsonic, su- 
personic, or alxed outflow boundary or a subsonic 
Inflow boundary. 

The physical space grid has the following prop- 
erties: one set of grid lines are straight and In 

the y direction with arbitrary spacing In the x di- 
rection; the second set of grid lines approzlastely 
follow Che wall and centerbody contoura; Che dy 
spacing of these grid lines Is arbitrary at one x 
location and cs proportional to those values at any 
other X location. 

The X, y physical space Is sapped Into a rec- 
tangular c, n coaputaclonal space as shown In n«. 

2. The sapping la carried out In two parts - Che 
first part aay j Che physical space to a rectangular 
coaputaclonal space while the second saps Che varia- 
ble grid coaputaclonal space to a uulfors grid coa- 
putatlonal space. Both the upper sud losser dual 
flow space walls collaose to the the saae grid line 
In Che coaputaclonal space, as shown In Fig. 2. The 
flow variables at Cue g.'ld points on the upper dual 
flow space wall are stored in the regular solution 
array while the vcrlables at the lower dual flow 
space wall are scored In a duaay array. These flow 
variables are coutlnually switched between chess two 
srrsya during the calculation. F^r details of the 
trsns forma dons, see Ref. 2. 

Numerical Method 

The coaputaclonal plane grid joints are divided 
up Into interior and boundary polr.cs. The boundary 
grid points are further divided up Into left bounda- 
ry, right boundary, wall, centerbody, wad dual flow 
space wall points (see Fig. 2). 

Interior Grid Points 

The Interior grid points are computed using the 
12 

unspllt MacCormack scheme. The governing equa- 
tions are left In nonconservaclon form. In order to 
taprove the coaputaclonal efficiency for high 
Reynolds number flows, Che grid points In fine part 
of the grid nay be subcycled. This Is accomplished 
by first computing the grid points In Che coarse 
part of the grid for one time step at. Next, the 
grid points In the fine grid are calculated k tlmea, 
where k Is an Integer, with a time step At/k. The 


grid points at the edge of the fine Rrld require a 
special procedure, because one of their neighboring 
points la calculated as part of the coarse grid. 
Zzeept for Che first subcycled tiae step, this point 
Is unknown. However, the values at C and t-t-At are 
known from the coarse grid solution and, so, the 
values between t and tw-At are determined by linear 
Interpolation. 

In order to further Improve the cooputatlona. 
efficiency, a special procedure Is employed to In- 
crease Che allowable clae step In the subcycled par* 
of the grid. This procedure allows the removal of 
Che sound speed from the clae step C-F-L condition. 
Procedures that accoopllahed this heve been proposed 
13 14 

by Harlow and Aasden and NacCormsek. The proce- 
dure of Harlow and Aasden Is an Implicit scheme that 
removes the sound speed. In both the x and y dlrec- 
tlona, by an Implicit treatment of the asas equation 
and the pressure gradient terms In the aosantum 
equations. NacCoraack'u procedure Is explicit and 
removes the sound speed In only one direction. 
(HacCoraack's procedure also Includes an Implicit 
procedure to remove the viscous dlffuslun restric- 
tion froa the time step C-F-L condition.) Because 
explicit schemes art easier to program for efficient 
computation on yxetor coaputers and because high 
Reynolds number flows usually required fine grid 
spacing In only one direction. It was decided to use 
a procedure alallar to the* of NacComack. 

HacCoraack's procedure Is based on the assuap- 
t'-on that the velocity component. In the coordinate 
dl. action with the fine grid spacing. Is negligible 
coapared to the sound speed. This allows the gov- 
erning equations to be simplified. NacCoraack then 
applies the Method of Characteristics to these sim- 
plified equations. However, for flowc over bodies 
with large amounts of curveture as well as many free 
sk^.r flows this assumption Is questl>*oable. Be- 
cause VNAP2 Is Intended to be s general code for 
solving a wide variety of problems, this assumption 
was felt to be too restrictive. Therefore, the aaln 
difference between HacCoraack's scheme sud the one 
presented below Is the removal of this restriction. 
This procedure, however, doe* assume chat rhe flow 
In Che y direction Is subsonic. 

The procedure here Is to separate Che governing 
equations Into two parts. The first part consists 
of the Mach line characteristic coapaclbility equa- 
tions while the streamline coapaclbility equations 
and all viscous terms aske up Che second part. 
Because Che sound speed limitation la due to Che 
first part, the second part can be computed by the 
standard MacCormack scheme. The first part of the 
governing equations Is 'olved by a special proce- 
dure. This special procedure consists of selecting 
an Increased time step based on Ay/lvl Instead of 
Ay/l v|-fa) where v Is the velocity component In Che 
y dc.'ecClor. and a Is the sound speed. The Mach line 
characteristics are then extended back froa the so- 
lution ptlnt at Che Increased time step. These 
charactei Istlcs will Intersect the previous solution 
plane outside Che computation star of chr MacCormack 
scheme. By calculating these characteristic Inter- 
sectlng points, special differences using this In- 
creased domain of dependence can be determined and 
used In Che first part of Che governing equations. 
These special differences use a larger Ay Chan Che 
MacCormack scheme and, therefore, are less accurate. 
However, thia procedure Is only used In boundary or 
free shear layers where the viscous terms dominate 
the solution. This special procedure decreases Che 
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coaputatlonal tlae for high Reynolds number flows by 
factors of $ to 10 over flows that are not subcy- 
cled. For details of this procedure, see Ref. 2. 

Left Boundary Grid Points 

The left boundary can only be an Inflow bounda- 
ry. For supersonic Inflow all flow yarlables are 
specified. For subsonic Inflow, there ar two dif- 
ferent boundary condition options. The flist speci- 
fies the total pressure pu, total temperature T_, 

* IS ” 

and fl«f angle 9 as proposed by Serra. The second 
condition speclfl' < the * and y ycloclty components 
u and V, respect rely, along with the density p and 
was shown to be correct for a wcll-poaed problem by 

Oliger and Surjstrom.^^ For a discussion of the 
merits of these t»^ boundery conditions, 
see Ref. 2. Following the Ideas of Norettl snd 

Abbett,^^ all the unap«:clfled deoendent variables 
are computed using a reference plane characteristic 
scheme. The trlscous terms are created as source 
terms. For mixed subsonic-supersonic Inflow, VNAP2 
checks Che Mach number at each grid point to deter- 
mine the corruct boundary condition. The u, v, and 
0 boundary coitdlclon includes s nonref lecclng option 
to eliminate the trapping of waves In subsonic, 
steady flows (see Ref. 2 for details). 

Right Boundary Grid Poi .cs 

The right boundary car be a supersov.lc outflow 
boundary or a subsonic lnfl>nr/outf low boundary. The 
subsonic Inflow option Is required for c/ises with 
fl<w separation at the right boundary. For super- 
sonic outflow, all the variables are extrapolated. 

For subsonic outflow, the static pressure p Is spec- 
ified and the remaining variables are calculated us- 
ing a reference plane characteristic scheme. If 
subsonic reverse flow occurs at the right boundary, 
Inflow boundary :ondltlons sist be specified. This 
Is accomplished jy leevlng p equal to the specified 
exit pressure an<f specifying p and v. This inflow 
boundary condition Is also discussed uy Oliger and 

Sundstrom.^^ These boundary conditions Include the 

nonreflecting procedure of Rudy and Strlkwerda. 

For mixed subsonic-supersonic outflow, VNAP2 checks 
the Hach re as her to determine the correct boundary 
condition. 

Wall Grid Points 

The wail boundary can be a free-sllp boundary, 

* free- Jet boundary, a no-sllp boundary, or a arbi- 
trary Inflow/outflow boundary. For the free-slip 
option, the wall slope Is the boundary condition and 
the remaining variables are calculated using a ref- 
erence plane characteristic scheme. For the free- 
Jet boundary option, the static pressure Is speci- 
fied and the code determines the free- Jet boundary. 
For the no-sllp bounda.y, the velocity components 
are set to zero while either the temperature Is 
specified or the temperature gradient Is set to zero 
(adiabatic wall). The density Is calculated by the 
*’*f**'*"ce plane characteristic scheme. For the ar- 
bitrary Inflow/outflow boundary, the static pressure 
Is specified. If the flow across the boundary Is 
outflow, the remaining variable are determined by 
the reference plane characteristic scheme. If In- 
flow occurs, the velocity component tangent to the 
boundary and the density are specified while the 
normal velocity component is determined using the 


reference plane characteristic scheme. A nonre- 
flecting boundary condition option Is Included. 

Centerbody Grid Points 

The centerbody boundary can be a free-sllp 
boundary, a no-sllp boundary, or a plane (axis) of 
•yametry. The free-sllp and no-sllp boundary calcu- 
lations follow the wall procedure. For flows where 
the centerbody Is a plane of symmetry, the grid 
points are computed by the Interior point scheme. 

The boundary condition la the requirement of flow 
symmetry. 

Dual Flow Space Wall Grid Points 

The dual flow space walls can be either a free- 
sllp or no-sllp boundary. The calculations follow 
the wall and centerbody procedures. 

Steady State Acceleration for Subsonic Flow 

because signals propagate In all directions In 
subsonic flows, disturbances can reflect around In- 
side the computational grid for aany time steps. 

reflection of disturbances can significantly 
prolong the convergence to steady state. Several 
^Iffrrent procedures for accelerating the conver- 
gence to steady state for both the py, T.^, 9 and u, 

V, p Inflow boundary condltlona are presented In 
Ref. 2. One technique that works well for very com- 
plex flows Is an extended Interval time smoothing 
ptocedurc. Here, the solution for all dependent 
variables on the first time step Is stored. The 
pressure at a specified grid point Is then monitored 
on each time step. When this pressure changes di- 
rection, the solution at the current time step Is 
F’t**'*8*d with the solution at the first time step. 
This averaged solution replaces the current time 
step solution and. In addition. Is stored In place 
®f the first time step solution. This averaging 
procedure Is continued until the flow Is steady. 

The results for subsonic, steady flow In a converg- 
ing duct are shown In Fig. 3. The top curve Is for 
a calculation in which the Initial-data surface con- 
sisted of stationary flow at the stagnation pressure 
and temperature. At time equal to zero, the outflow 
pressure was dropped from tne stagnation value to 
the desired value, thus simulating a bursting dia- 
phragm. The middle curve Is for a calculation In 
which the Initial-data surface was the 1-D solution 
generated by the VNAP2 code. The bottom curve shows 
the calculation emplovlng the l-O Initial-data sur- 
face and the extended Interval time smoothing. All 
three aolutlona employed the p.j., T.^., and 9 Inflow 

boundary condition. From Fig. 3, we ace that both 
Improving the accuracy of the Initial-data surface 
and employing the extended Interval tine smoothing 
significantly Improved the convergence to a steady 
state. As a result, both procedures are utilized In 
the following results. 

Results 

The results presented here are for six high 
Reynolds number flows; one Internal, two external 
and three Internal /external cases. One case In each 
category has flow separation and the last case in- 
cludes transition rrom laminar to turbulent flow. 
These cases represent very complex flows and as such 
were selected to Illustrate the more difficult prob- 

that VNAP2 Is capable of solving. There Is s 
considerable number of less complex viscous as well 
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•s Invltcld flows chan VNAP2 can solve aore accu- 
rately and with significantly Less aaounts of coapu- 
cer tiae All cases were run rich the unaoiifled 
VNAP2 code utilizing orl, • »aall data file. At 
this tlae, very few psraaetrlc studies to determine 
Che optiauB turbulence aodel paraaecers, initial- 
data surface quantities and grid point distribu- 
tions, have been carried out* As a result, the ac- 
curacy and efficiency of these results do not neces- 
sarily represent the optlaal use of the VNAP2 code. 

Internal Flow 

The internal flow case is nozzle B-3 of Ref. 19 
and is Che planar, converging-diverging nozzle shown 
in Fig. 4. The flow is froa left to right with the 
physical apace grid enclnned by the dashed line. 

The Reynolds nuaber based on the throat height is 

7.7 > in^. At the left boundary, p^ is set equal r.o 

200.6 kPa (29.1 psla), T.^ is set equai to 300 R and 

a is sec equal to 0. At the right boundary, extra- 
polation is used when the flow is supersonic. When 
the flow H subsonic outflow, p in set equal to 
101.4 kPa (14.7 psia). For subsonic inflow, in ad- 
dition to specifying p, v is set equal to 0 and p is 
set equal to Che average of Che Hall and aidplane 
values. The wall is a no-sllp boundary and the re- 
sults presented liere employed the cw > equation tur- 
bulence The Inlclsl condlcitna consisted of 

l-D, Invlscid flow that was sonic at the throat end 
aubsonlc downstream. 

The physical space grid, Mach number and turlni- 
lence energy contours .re shown in Fig. 5. The ex- 
per'*aental data are from Ref. Iv. From Fig. S, we 
see that at this pressure ratio, the flow separated 
creating a reverse flow region. The wall and mid- 
plane pressures are shown in Fig. 6. This calcula- 
tion used a 45 by 21 grid and required 3(HX) time 
steps f30,000 subcycled time steps) and 2.1 hours of 
epu time (CDC-7600) to reach steady state. 

rxtemal Flow 

The external flow cases are conflgs. 1 and 3 of 
Ref. 20 and are Che axlsyMecrlc, boaccall afterbody 
flows, with solid bodies simulating the exhaust )et, 
shown in Fig. 7. Tie flow is from left to right 
with the physical space grid enclossd by Che dashed 
line. The Reynolds number baaed on x at Che laft 

boundary is 1.05 x 10^. The first case (1 ■ 27.0 
cm) is config. 3 while Che second (t ~ 12.2 cm) is 
ccnflg. 1. Both cases consisted of Che some flow 
conditions. The left boundarv Inflow profiles of p.^, 

and T^ for a free screao Mach nuaner of 0.8 were de- 
termined using the same invlscld/boundary layer pro- 
cedure em"loyed by Ref. 21. The flow angle 3 wts 
set equal to 0. At Che right boundary p was set 
equal to the free scream value. The wall is an ar- 
bitrary inf low/oucf low boundary. For outflow, p is 
set equal to Che free stream value. When inflow oc- 
cu'‘s, in addition to specifying p, u and p are sec 
equal to their free stream values. The centerbody 
is a co-slip boundary. Both '.rulatlons employed 
Che extended Interval time smooching. The initial 
conditions consisted of extending Che Inflow pro- 
files downstream to Che right boundary. 

The physical space grid, pressure and Mach num- 
ber contours for config. 3 (t " 27.0 cm), employing 
the mixing-length turbulence s»del, are shown in 


Fig. d while Cl e surface pressure is shown in flR- 
9. r..e physical space grid, Mach number and turbu- 
lence energy contours for config. I (t ■ 12.2 cm), 
employing Che two equation Curbulencr aodel, are 
shot.n in Fig. 10 while the aurface pressure is shown 
in Pii:. 11. The experimental data, for both cases, 
are fiia Ref. 20. Proa these figures, we see that 
the flow for config. 3 remained attached while sspa- 
ratlon occured for config. 1. The mixing-length 
model produced slightly better results for config. 

3. However, the two eq>iation aodel more accurst* ly 
predicted the pressure plateau in config. 1, but un- 
derpredlcted the amount of upstream expansion. 

Swanson^^ found that s relaxation or lag model im- 
proved Che pressure plateau prediction of the mix- 
ing-length aodel, but at the expense of also under- 
predicting Che amount of upstream expansion. The 
config. 3 calculation, employing Che mixing-length 
aodel, used a 40 by 25 grid and required 750 tlae 
steps (15,OOC subcycled time steps) and 1.0 hours of 
CPU time (CDC-7600) to reach steady state. The 
co.nflg. 1 calculation, employing the mixing-length 
model, used u 47 by 29 grid and reoulred 750 time 
steps (17,000 subcycled time steps) and 2.1 hours of 
CPU time (CDC-7600) Co reach steady state. The two 
equation aodel computational Claes were 1.4 hours 
for config. 3 and 3.7 hours for config. 1. 

Internal /External Plow 

The first two cases are Che tvo external flow 
cases presented above, hut with Che solid simulators 
replaced by the exheust Jets. The geometry Is shown 
In Pig. 12 with Che physical space- grid enclosed by 
Che dashed lines. The left external boundary, wmli. 
and rlehc boundary are the same as Che external 
cases. At Che left internal boundary, Is sec 

equal to 132.4 kPa (19.2 psla), T^ Is set equal cr 

300 K and 6 Is set equal Co 0. The free stream 
pressure ^s 65.2 kPa (9.45 psla). The centerbody is 
the flow centerline, while the dual flow space malls 
are no-sllp boundaries. Again, both calculatlo .s 
employed the extended interval time smoothing. Only 
the two equation turbulence model was employed for 
these t»x> cases. The initial conditions for these 
cases consisted of the external flow solutions pre- 
sented above along with the l-D, invlscid flow solu- 
tion for the nozzle. 

The physical space grid and Mach number con- 
tours for confiq. 3 (t ^ 27.0 cm) are shown in Pig- 
13. The external surface pressure is shown in Pig. 
14 while Che total pressure profiles for the shear 
layer, produced by the interaction between Che ex- 
haust Jet und the external flow, are shown in Fig. 

15. The physical space grid and Mach number con- 
tours for config. 1 (t ■ 12.2 cm) are shown In Pig. 

16. The external surface pressure is shown in Pig. 
17 while Che total pressure profiles for Che shear 
layer are shown In Pig. 18. The experimental data, 
for both cases, are trom Refs. 20 and 22. Proa 
Pigs. 14 and 17, we see Chat Che computed solutions, 
for both cases, undsrpredicted the shear layer 
spreading race. The same trends were fo*'nd in re- 

23 

suits generated by t patched method. The config. 3 
calculation used a 40 by 38 grid and required 300 
time steps (40,000 subcycled time steps) and 5.6 
hours of CPU time Co reach steady state. The con- 
fig. 1 calculation uaed a 47 by 42 grid and required 
300 time steps (40,000 subcyclcd time steps) and 
9.C hours of CPU time (CDC-7600) to reach steady 
state. 
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The l^et caee Is the NACA 1-64-100 Inlet of 
Ref. 24 and Is .'hown In Pig. 19. ine flow Is froa 
left to right with the physical space grid enclosed 
by the dashed line. The Reynolds nuaber based on 

the aaxlaia external dlaaeter is 6.1 « 10^ . At the 
left boundary, was set equal to 101.4 kPa (14.7 

psla), Xj was set equal to 294.4 K and 6 was set 

equal to 0. The free atreaa pressure Is 66.5 kPa 
(9.64 psla) which produces a free screaa Hach nuaber 
of 0.8. The wall is the arbitrary inflow/outflow 
boundary and uses the saae boundary conditions as 
the crevlous two cases. At the right external 
boundary, the pressure was set equal to the saasured 
values. In the experlacnts of Ref. 24, the Internal 
flow rate was controll'.d by a throttling aechanlsa 
well downstreaa of the right boundary. In order not 
to coapute this rtther extensive flow region, s 
value of pressure was specified at the right inter- 
nal boundary sue^ that the fnvlscld, 1-D aaas flow 
equalled the experlaental value. ‘Hie centerbody Is 
the flow centerline, while tlie dual flow space walls 
are no-sllp ^ idarles. The extended Interval tlae 
saoothlng was eaployed. The Initial conditions con- 
sisted of 1-D, Invlscld flow. 

The physical apace grid, Mach nuaber and turbu- 
lence energy contours are shown In Pig. 20, while 
the surface pressure Is shown In Pig. 21. The ex- 
perlaental data are from Ref. 24. Proa the turbu- 
lence energy contours In Pig. 20, we see that the 
flow over the Inlet Is Initially laalnar, but quick- 
ly becoaes turbulent. The flow transitions first on 
the Interior surface. Reference 24 did not give the 
transition locations. The surface pressures In Pig. 
21, are for the first 62 of the Inlet. This saell 
area of Interest, coablned with the thin laalnar 
boundary layer, required s very fine grid spacing In 
both coordinate directions. Because of the fine 
grid spacing In the x direction, the subcyllng op- 
tion was not used. The differences in surface pres- 
sure, between theory and sxperlaant, at the Inlet 
tip are probably the result of too coarse s grid. 

The difference, between theory and experlaent, near 
the Internal, right boundary Is aoat likely due to 
the approxlaste treataent of Che Internal flow at 
this boundary. This calculation was slao oade using 
Che alxlng- length turbulence aodel, however, this 
aodel significantly overpredlcts Che level of turbu- 
lence upstreaa of Che Inlet. This Is because the 
downstreaa blockage due to the presence of the Inlet 
creates a weak shear layer profile upstreaa of the 
Inlet Chat extends a large distance in the cross 
screaa direction. Therefore, the alxlng length aod- 
el predlcCb very large alxlng lengths. This pro- 
duces turbulent viscosities, upstreaa of the Inlet, 
chat are several orders uf oagnltude larger chan the 
aolecular value. As a result, Che alxlng-length aodel 
solution Is not presented here. This calculation 
used a 53 by 44 grid and required 3000 tlae steps 
and l.O hours of CPU tlae (CDC-7600) Co reach steady 
state. 

Conclusion s 

A genersl, user oriented coaputer prograa for 
coapuclng high Reynolds nuaber flows has been pre- 
sented. Six high Reynolds nuaber flow calculations 
were described. These coapuced results show that 
practical Navler-Stokes applications are possible. 
However, these results also Indicate the need for 
better turbulence aodellng for separated flows and 


aore efficient solution algorlthas for verv nonunl- 
fom grid point distributions. 
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Fig. 1. Typical Internal, external and Intemal/cxtemal geometries. 


6 



Wall 




UpfMT Dual Row Spoca Wall 


Lowor Dual Row Spoca Wall 


Contorbody 


U; 


Computational Spoce 

Fig. 2. Physical and coopucational spacv's. 
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Fig. S. Internal flow physical space grid, 
Mach number and turbulence energv 
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Fig. 3. Pressure ratio vs time step for 

subsonic flow in a converging duct. 
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Fig. 7. External flow geometry 
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Fig. 8. Extemtl flow (config. 3) physical 

space grid, pressure and Mach number 
contours . 
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Fig. 10. External flow (conflg. 1) phyalcal 

apace grid Mach number and turbulence 
energy contours. 


Fig. 13. Internal /external flow (conflg. 3) 

physical space grid and Mach number. 
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